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bation theory. The process is dominated by the 0(j> 6 ) in the momentum 
expansion where tree-level amplitudes from the effective Lagrangian £ 6 en- 
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C^~- , from C 2 . We estimate the 6 independent £ 6 coupling constants by resonance 
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. of the uncertainties involved. 
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I. INTRODUCTION 



The reactions involving two neutral pseudoscalar mesons and two photons have received 
much attention during the past few years. In the framework of chiral perturbation theory, 
they are unique in that the tree-level amplitudes from the Lagrangians C 2 and £ 4 both 
vanish. As a consequence, processes such as 77 — > tt tc , 77 — > 7r°77 (or the reactions 
related to them by time reversal) test higher order and loop contributions in the chiral 
perturbation series. In fact it turns out that the 0(p 4 ) one-loop amplitude for 77 — > 7r°7r° 
]T|fJ] underpredicts the Crystal Ball data || at low invariant mass of the two photons (where 
chiral perturbation theory is supposed to work best). The 0(p 4 ) result for the i] — > 7r°77 
decay-width is a tiny 3.89 • 10~ 3 eV [|[], to be compared with an experimental value of 
0.84 ± 0.2 eV |J. The //-decay amplitude is thus dominated by higher order contributions. 

According to Weinberg's power-counting formula M, the (9(p 6 )-amplitude will contain 
tree-level contributions from the effective Lagrangian C 6 together with one-loop contribu- 
tions from £ 4 and two-loop contributions from C 2 . The form of the complete £ 6 -Lagrangian 
has been given recently ]7J ; we will see below that a total of 6 linearly independent structures 
contribute to the processes considered here. In extending the chiral Lagrangian approach 
to 0(p 6 ), we therefore have to deal with two difficulties: A realistic estimate of the 6 low- 
energy constants involved and the solution of the 2-loop Feynman integrals, in particular 
for the case where the two loop momenta overlap. 

There have been various attempts to determine the £ 6 7/-decay amplitude in the past 
either phenomenologically or in a bosonized Nambu-Jona-Lasinio model. In [|J, only the 
leading terms in a l/iV c -expansion are considered, i. e. the Lagrangian is restricted to the 
single flavour trace terms 

d 1 F^F^Tr(Q 2 5 x U5 x U^) + d 2 F m F^Tr(Q 2 5 a U5^). (1.1) 

(For the notation, see Section [II A|) . If the two coupling constants d\ and d 2 are determined 
by comparison to a vector dominance (VMD) amplitude 

-d x = irf 2 ~ 3.6 ■ 10" 3 GeV" 2 , 

the tree level amplitude from the 0(p 6 ) Lagrangian ( |1.1| ) yields a total decay width of 

r 6 (r? -> 7T°77) = 0.18 eV. 

This is to be compared to the full VMD amplitude of flJ^JlOU 

r WD (^/ T 7) ~0.31eV. 

In the Nambu-Jona-Lasinio model, a further single trace term contributes to £ 6 , namely 
TQM 



d,F, u F^Tr{Q 2 { X ] U + U^x))- (1-2) 



The determination of the constant d^ is ambiguous. According to Ref. ||11|| , the term ( |1.2| 
increases the f?-decay width to 



^ENJL 



(r] 7r°77) = 0.58 ± 0.3 eV 
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whereas the authors of |12| advocate values of 

r NJL (r) -»• tt°77) = 0.11... 0.45 eV. 



If d 3 is calculated via resonance saturation in 77 — > 7i°7i° JT3J, the term ( |1.2j ) doesn't affect 
the ?7-decay width very much. Conversely, an attempt to keep d\ and d>2 as given by vector 
meson resonance saturation and then fit d% to the experimental value of the 77-decay width 



14| leads to a very large d% and implies that further terms in the Lagrangian are required for 
a realistic description of pion polarizabilities and the the process 77 — > 7r°7r°. In summary, 
one is led to the conclusion |L^] that a self-consistent, quantitative description of the decay 
i] — > 7r°77 based on the Lagrangians (|1.1| ) and (|1.2| ) is problematic. 

The numerical decay widths quoted so far are calculated from pure £ 6 tree-level am- 
plitudes. A complete treatment of the £ 2 two-loop contributions has been achieved in an 
SU(2)-calculation of 77 — > 7r°7r° |]TB[, where only pionic degrees of freedom are taken into 
account. Numerically, the two-loop amplitude is crucial for a fit of the experimental spec- 
trum. For the 77-decay, the contributions of the factorizable 2-loop diagrams have been 



given recently |TI|. Although the effects on the 77-decay are small, it is clearly desirable to 



carry out a complete calculation of the two neutral meson processes at two-loop level. 

The purpose of the present paper is thus twofold: First, we explore the structure of 
the tree-level amplitudes if the most general Lagrangian C§ is used and compare the result 
to previous calculations. Second, we outline a method to derive the two-loop amplitude 
including overlapping graphs. We present the results for the case of pure pion loops and 
estimate the missing parts for 77 — > 7r°77. A method of solving the specific Feynman inte- 
grals occurring in the problem may be of general interest and is therefore discussed in the 
appendix. 



II. THE 77 7T°77 AMPLITUDE 



A. General Formalism and Kinematics 

Throughout the following, we will use the notation of applied to our special case: The 
[/( 3) -matrix U = exp(-^($g + $1)) contains the nonet of pseudoscalar mesons 



$ 8 (x) 



V V2K- 



v / 2tt + V2K + \ 
-tt° + ^ 
V2K° 



-2r? 8 , 
73 / 



*,.(,•)= \/^/,/ (2.1) 
The covariant derivative for photons as 



F n = 93.2 MeV is the pion decay constant | 
external gauge fields reads 

D„B = 6„B + ieA^Q, B] 

where Q = diag(2/3, — 1/3, — 1/3) is the quark charge matrix and B stands for be any 
operator that transforms linearly under the chiral group U(3)l x U(3)r: B — » B' = VrBV^ 
We note that 
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vr+ K+ 

[Q,$]= | 7T- 1 (2.2) 

ir~ 

so that the commutator [Q, $] couples the photon to charged mesons only. In Eq. flLlQ , 
coupling of the neutral mesons to photons is mediated by the field tensor F^ = {S^Ay — 

The physical rj particle can be written as || 

2 r 1 

rj = cos % - sin 9r)i ~ - V2r] 8 + -rfr. (2.3) 
In our normalization, the lowest order chiral Lagrangian reads 

C 2 = ^Tr(D„UD^) + ^fTr(x^U + tf X ) (2.4) 

where the mass terms contain the quark mass matrix \ = 2-Bodiag(m M , m^, m s ) and the 
constant B relates the masses of quarks and pseudoscalar mesons via 

B = < = rnl = Am 2 K (25) 
m u + m s m u + m d m d - m u 

A numerical estimate for the electromagnetic mass split of the kaon yields the value of 
Am K ~ 6000 MeV 2 §. 

The general form of the amplitude for the decay rj(P) — > Tc°(p)j(qi, £1)7(92, £2) is jlOlJI^ 

+B(s, t)- g ■ (g» v P ■ q x P ■ q 2 + P*P v t _ $p> p . Ql _ p^p . g A | ^ ( 2 . 6 ) 

where we have defined the kinematic invariants s = {qi + q2) 2 i t = (P — q^) 2 , u = (P — qi) 2 = 
m 2 + m\ — s — t. The decay rate is calculated according to 

fs = £ * (l AS - < B I 2 +^('»X - tu?) , (27) 

where s is restricted to < s < (m 2 — m^) 2 and 

1 



^1 2 — — 
M 2 



m^ + m,- s) ± J (m 2 + m 2 - s) 2 - 4m 2 m 2 . (2.8) 



B. Complete £ 6 Amplitude 

The tree-level amplitude from the Lagrangian £ 6 forms the dominant contribution of 
the 77-decay width and will therefore be discussed separately from the one- and two-loop 
corrections. 
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The most general chiral Lagrangian at 0(p) 6 has been derived for the SU (3)-case 
U = exp(-^$g). Introduction of 77 — 77' mixing according to Eq. ( |2.1|) eliminates a symmetry 
constraint and potentially generates more terms that we have to neglect at this point (for an 
example see ]16)). Keeping this simplification in mind, we can take the £ 6 -terms required 



from Table II, section no. 4 (terms involving B 3 q — B50) of Ref. 0j. The number of indepen- 
dent structures contributing to the 2 photon/2 neutral meson processes can be determined 
in two steps: First make use of the simplifications explained in Sections IV B/C of Ref. J7[ 
and note that, owing to Eq. (|2.2| ), the commutator [Q, $] can be set to zero in the covariant 
derivative as well as under the trace of the £ 6 -expressions. As a result, for instance both 
the terms in B 35 and B 3e of Ref. J7|, Table II reduce to the G?2~term in Eq. ( |1 . 1| ) : 



B 3B Tr([D a U]4Df > U].[G^] + [G^]+) + B 36 Tr([D a U]4D^U]4G^] + [G 



+ / 



d 2 F lMa F^Tr(Q 2 5 a U5 p W), d 2 = 4(535 + B m ). (2.9) 

The terms resulting from this procedure are more symmetric than the original general £ 6 - 
terms. This implies that terms which are independent in the general Lagrangian C 6 can be 
redundant due to [Q, $] = 0. Therefore, in a second step, we have to check the result of 
the simplified Lagrangian and reject terms that are linked by a trace relation of the type 
discussed in ]7), App. A. For example, one of the terms 

B u Tr([D^UUD^UUG^} + [G af3 ] + ) = AB u F aP F^Tr{Q 2 5^U5^), 
B^Tr{[D^U\.[D> i U\J)Tr{[G^\ + [G ot p] + ) = 4B 37 F aP F al3 Tr(Q 2 )Tr^U6^), 
B^Tr{[DJJ}-[G a P} + )Tr{[D»U\-[G al3 } + ) = 4B 38 F a0 F^Tr(Q5^U)Tr(Q5^) (2.10) 

is redundant after simplification due to the the trace relation 

ATr(A 2 B 2 ) + 2Tr(ABAB) - Tr(A 2 )Tr(B 2 ) - 2(Tr(AB)) 2 = 

for any complex 3 x 3-matrices A, B with Tr(A) = Tr(B) = 0. Alternatively, the -B34 and 
B 37 can be seen to be equivalent because of Tr(Q 2 ) =2/3. 

We are finally left with the 3 single trace terms of Eqs. (|1 . 1|) and (|1.2j) plus 3 double 
trace terms 

dtF^F^TriQ^USxtf) + d 4 F^F^Tr(Q8»U)Tr(Q5^) 
+ d 2 F llQ F^Tr(Q 2 5 a U5^) + d^ a F^Tr{Q5 a U)Tr{Q5^) 
+ d 3 F, u F^Tr(Q 2 ( X j U + tf X )) + d,F^Tr{Q 2 )Tr{ X ] U + U^ X )). (2.11) 

and thus a total of 6 low-energy constants to be determined either from experiment or in 
some modelQ. The low-energy constants d\, . . . , d§ are related to the S-constants of 0] as 
follows: 



Elimination of redundant terms in the Lagrangian is desirable for simplicity but not indispensable 
from a practical point of view. Keeping more terms would result in fixed linear combinations 
of coefficients appearing in the amplitudes; the fit procedure would ultimately not be affected. 
Moreover, we have some freedom as how to represent the Lagrangian and which terms to discard. 
A choice different from Eq. ( |2.11| ) might affect the expressions for the amplitudes and the values of 
the parameters quoted but of course must leave the physical content of the Lagrangian invariant. 



5 



d x = AB 34 +3^7, d 4 = AB 38 , 
g 

d 2 = 4(5 35 + B m ) + -B 39 , d 5 = 4(S 40 + B 41 ), 

d 3 = 4B 47 , ' d 6 = 4B 50 . (2.12) 

Using Eq. (|2.3| ), the 77-decay amplitude from the full £ 6 -Lagrangian can now be cast in 
the form (12.0) with 



9 /9 

A v (s,t) = — ^ {(4d 1 - Yld 4 )m\ - 2d 3 ml ~ ( 4rf i ~ 12rf 4 + d 2 - 3d 5 )P ■ {q x + q 2 )} 
4\/2 s 

B v (s,t) = ^L^-{d 2 -3d 5 } (2.13) 
Note that we used the isospin approximation m u = neglecting a small contribution from 



the double trace mass term in Eq. (|2.11| ) proportional to d§. This does not change any of the 
conclusions of this chapter. The ofa-ternis does affect, however, the 77 — > 7r°7r° amplitude 
which can be written in analogy to Eq. (|2.13|) 

4 5 1 5 4 5 12 

A n (s,t) = - — (s-2ml)(-d 1 -d 4 ) - — s (-d 2 -d 5 ) - ^(g4+y4) 

B^t) = ^- 2 S -{\d 2 -d,). (2.14) 

Determination of the 6 constants d\ . . . d e from experimental data alone is now impeded by 
a combination of factors: 

1) For the process rj — > 7r°77 where the £ 6 -amplitude ( p.!3|) forms the leading contribu- 
tion, only the total decay width has been measured |J. 

2) The energy dependence of the cross section for neutral pion pair production is experi- 
mentally known. In this case, however, the amplitude ( [2.14j ) is small and interferes with the 
large, complex one-loop and two-loop contributions |13j so that the fit becomes ambiguous. 
Complementary data such as pion polarizabilities [I3|,|r7] don't put sufficient constraints on 
the parameters. 

3) There are further two neutral meson/two photon data, e. g. the production cross 
sections for 77 — > r/n and 7 — > K°K° J|, but the processes have too high energy thresholds. 

On the other hand, the method of resonance saturation has been successfully applied 
in order to confirm the phenomenological constants of the Gasser-Leutwyler Lagrangian £ 4 
18| . Moreover, the pion production amplitude ( |2.14|) with the constraint d 4 = c?5 = d% = 



has been determined via resonance saturation [|T3[ and compares well with the data. We 
proceed to discuss how the method can be applied to our general case. 



C. Meson Resonance Amplitude 

In the VMD model, the leading contribution to i] —>■ tt°77 is generated by exchange of 
internal vector mesons u, p. This generates the amplitudes ||L0|,|4l,|l3[1 : 
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< MD (M) 



G 



B 



VMD, 



S,t) 



E 

V=w, 



E G ^ v 2 

V=u),p ~~ 



t + mi ^ u + mj 



t — IJly 



t — m v 



u — m v 
1 



u 



m\ 



E 



G 



V=u>,p,<5> m V 



E 

V=LU,p 



Grjy ■ S 

m v 



(2.15) 



The right hand side indicates the 0(p 6 ) low energy limits of the VMD amplitudes. The 



VMD amplitude for neutral pion pair production 77 — > 7r 7r° is (|2.15 ) with the substitutions 



G v v — > G-xv, m ri m ir- The coupling constants G v v, G n y, can be extracted from the 
decay widths of the vector mesons. We list the experimental values |9.13] of all the coupling 
constants used in the normalization of Eqs. ( p,15| )- (|2.18| ) together with the meson masses 
in Table |. 



The contributions of the C-odd axial-vector resonances are 19 



B*(s,t) 



E 



G 



T),B 



—t + mi 



G v v- 



B=b u h,h' 



t — m 
1 



+ 



+ 



—u + m. 



21 



u — m 
1 



_ ( 3m v + m l- s 

B=b u h,h' m v 



E 



t — m\ u — m 



v 



E 

B=bi,h,h' 



G n v ■ S 

my 



(2.16) 



Again, there is an equivalent amplitude for pion pair production. The axial vector mesons 
interfere constructively with the vector meson amplitudes (|2.15|) and thus enhance the 77- 



decay width while deteriorating the reproduction of the pion pair production data [|1J . We 
include only the measured b\ resonance in our model. 

The production cross section 77 — ► r/n is dominated by the (scalar) resonance a (983) 
and the (tensor) resonance a 2 (1318). Their contributions can be written as 



A T Js,t) 



G^ 



mi 



2 ■ 
mf — s 



G„t s 



for the 02 and 



% d (s + m 2 n - 3m 2 w ) + Gl m 2m\ 



Bi 



nit 



[s,t) 



2 m% — s 



0; G S J ■ Gt m > 



(2-17) 



Sm 



(2.18) 



for the a respectively |X8|J5|. The 77 — > 7r 7r spectrum shows the scalar resonance /o(983) 
and the tensor resonance ^(1275) with contributions analogous to Eq. Q2.17j ), (|2.18| ) (set 



m r 



and choose the appropriate couplings, see 



The coupling constants for the vector and axial-vector mesons are positive, but the signs 
of the (s-channel and therefore not quadratic) coupling constants G s and G T relative to 
the VMD amplitude are ambiguous. Our choice is motivated by the following observation: 
Consider the one-loop and (approximate) analytic two-loop amplitude for 77 — > 7r°7r° from 
Ref. [TI| as the Q CD-background for the process and complement this amplitude by the full 
resonance amplitude of Eqs. Q2.15|) , (|2.16|) , ( |2.17|) and (|2.18|) . At high energy, parametrize 
the resonance widths m? and m$ by a relativistic Breit-Wigner form given by Eqs. (5) and 
(6) of the original data analysis of the Crystal Ball experiment [ 20| . The result is of course 
not accurate in the region of the fo and f% resonances but fits the shape of the data fairly well 
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provided the signs of the couplings are chosen as in Table I (see Fig. ||). With the opposite 
signs, the measured cross section can not be reproduced. Our sign convention implies that 
the 77-decay width is increased by the contribution of the scalar and tensor mesons. By the 
way, this result is consistent with the interference pattern found for the scalar resonances in 



the extended Nambu-Jona-Lasinio model of 11 



Having constructed a realistic meson exchange model for the two neutral meson/two 
photon processes (see Fig. [I]), we can now proceed to test the prediction of the chiral 
Lagrangian C 6 to both pion pair production and the 77-decay width. To this end, the 
propagators of Eqs. ( 2.15| ) - ( |2.18|) are taken to first order in the kinematic invariants s, u 



and t as indicated in Eqs. ( |2.15| ) and ( |2.16| ); the resulting 0(p 6 ) low-energy expressions are 



used to fit simultaneously the parameters of the chiral amplitudes Eqs. (|2.13|) and Q2.14|) . 
The procedure yields a total 77-decay width of 

T L6 (77 -> vr°77) = 0.652 eV. (2. 19) 

(The last digit has been quoted for the discussion of the small loop contributions, see below). 
As shown in Fig. |2|, the spectrum obtained with the (9(p 6 )-fit differs considerably from the 
original meson exchange spectrum, being smaller as s — > and showing a broad maximum 
at a/s = _E 77 ~ I \JSmax . The difference at low invariant mass \fs is induced by the higher 



order terms in the t-channel amplitudes Eqs. Q2.15j ) and (|2 .16 ). As s grows, we increasingly 



neglect contributions from the s-channel amplitudes ( [2.17D and ( |2.18|) as well. It is in fact 
easy to perform an "all order" fit of the meson-exchange amplitude, see Fig. ||. We discard 
this option, however, as not consistent with the spirit of the momentum expansion. 

The 0(p 6 ) tree-level contributions are small for low-energy pion pair production so that 
details of the fit have only a minor effect on the cross section (see Fig. ^). 77 — > 7r°7r° 
therefore doesn't put further restrictions on the parameters of the model. Note, however, 
that it is in general not possible to obtain an accurate reproduction of the (all order) meson 
exchange amplitude by £ 6 -terms in the energy domain E yy > 500 MeV. 

The coupling constants di . . . d§ derived from the fit and collected in Table deserve some 
comments. At first glance, a comparison of the 0(p 6 ) and all order results for the kinetic 
constants d±, d 2 , d A and d 5 suggests that the fit procedure is highly unstable, changing the 
order of magnitude and even the signs of the constants as one passes from one scheme 
to an other. This behavior is indicative for the fact that the ratio between the empirical 
meson coupling constants G v and G n used in the fit (Table |) is close to the ratio between 
the SU(3)-factors contained in Eqs. fl2.13p and ( |2.14| ) - quark symmetry and therefore the 



physics of the process itself seems to exclude a stable inversion procedure. The results appear 
more trustworthy if one keeps in mind that in both amplitudes, only linear combinations of 
the coefficients d\ . . . d§ enter. For the 7/-decay (Eq. (|2.13|) ), these are 

a) d x - 3d 4 = 4.834 + 1^37 - 12^38 = -3.56 ■ 10~ 3 GeV~ 2 (-4.81 • 10~ 3 GeV~ 2 ) for the 
0(p 6 ) (all order) fit, to be compared with the value of —3.6 ■ 10~ 3 GeV~ 2 derived in the 
simplest fit scheme (Eq. 

b) d 2 - 3rf 5 = 4(5 35 + 5 36 ) + |5 39 - 12(5 40 + £41) = 9.59 • lO" 3 GeV" 2 (12.4 • lO" 3 
GeV -2 ) for the 0(p 6 ) (all order) fit, to be compared with the value of 7.2 • 10 -3 GeV -2 from 
Eq. |). 

The linear combinations |c?i — c? 4 and §g?2 — d 5 that appear the pion production amplitude 
are similarly stable with respect to the fit method chosen. Due to the similarity of our 
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meson exchange model with the one used in []13|1 , they are doomed to be consistent with the 
constants found there. 

The mass term constants d% and d§ show a slightly different behavior. As d§ doesn't 
influence the 7/-decay amplitude, they have to be stable separately, which is indeed what we 
observe. The values found for d% are consistent with the value fitted to the experimental 
total decay width Q (rf 3 = 45 ■ 1CT 3 GeV~ 2 ) within ~ 10%, furthermore, @ obtains a 
spectrum similar to the one shown in Fig. ^. Note, however, that we have derived our 
coupling constants from the general form of the Lagrangian C 6 and the meson exchange 
model alone: no assumption concerning the f?-decay width or pion production cross section 
was made. The agreement between our model and the one discussed in might indicate 
that the relatively high 7/-decay width measured is indeed correct. 

The mass term difference A m = |d 3 + f d 6 = l(5B 47 + 12B 50 ) = 1.44 • 1(T 3 GeV" 2 
entering the pion production amplitude ( |2.14|) is more than an order of magnitude below 
the value of d% (and thus of the order of magnitude of the kinetic constants). Again, our 
values of A m have to be consistent with the fit of |L3| by construction. However, identifying 
A m with the mass term d 3 entering the ^-decay amplitude would inevitably lead to an 
underprediction of the 77-decay width. 

It seems therefore that any attempt to describe 77 — > 7r°77 and 77 — > 7r°7r° simultaneously 
in a model containing only 3 £ 6 constants as in Eqs. (|1.1| ) and (|1.2|) is bound to be troubled 
by the incompatibility of the mass terms. This explains some of the inconsistencies about 
the constants d%, c?2 and <i 3 found in the literature. 

Because of to the relative insensitivity of the 77 — > 7r°7r° cross section, a measurement of 
the 77-decay spectrum would discriminate between different sets of £ 6 coupling constants. 
However, in order to pin things further down, we need an accurate estimate of the loop 
contributions generated in chiral perturbation theory. 



D. Effect of Chiral Loops to 0(p 6 ) 

The following section is devoted to a discussion of the size of the various loop contribu- 
tions occurring up to 0(p 6 ) and the basic technicalities of their evaluation. For the formal 
aspects of an two-loop calculation and the implications for the amplitude, we refer to earlier 
publications |T3| , |26| . 

The leading loop contributions appearing in chiral perturbation theory are the 0(p 4 ) 
one-loop diagrams with vertices generated by the Lagrangian C 2 of Eq. ( |2.4| ). A detailed 
calculation of this contribution for the 77 — > 7r°7r° amplitude has been given in |l]; the 
?7-decay can be treated analogously, employing expressions similar to Eq. ( |4.8| ) without the 
tadpole term. For the result see Ref. 00. The (9(p 4 )-amplitude involves charged pion and 
kaon loops and strongly interferes with the tree-level amplitude from £ 6 so that its effect 



2 From the r/7r°7r + 7r _ -vertex displayed below (Eq. ( [4.1|) ), we obtain only the first (leading) term 
of the one pion loop amplitude Eq. (14) of Ref. fl. We agree, however, with the numerical result 
Eq. (16) of this reference. 
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increases with the latter amplitude. Adding it to the £ 6 result Eq. ( |2.19|) , the decay width 
increases by about 10% to 

r i 6 +o(p 4 ) ^ = Q 733 eV> (2>20) 

We note that about 3/4 of the one-loop effect stems from charged kaon loops. Inclusion of 
pion loops alone would increase the decay width to 0.668 eV. 

At order 0(p 6 ), there are further one-loop diagrams generated by the Gasser-Leutwyler 
Lagrangian C A ||. Besides the recalculation of vertices for all the £ 4 -terms, these contri- 
butions introduce no new technical difficulties. For a pion loop only, the task is simplified 
because, like in the case of the Lagrangian C 2 , r]ii 7i + tt~ -vertices can only be generated 
by the mass terms. This leaves the terms in L4 and L5 as possible candidates where the 
renormalized coupling constant L\ is ^-suppressed and empirically consistent with 0. The 
/^-amplitude invokes the loop function F defined in the Appendix and is compact enough 
to be displayed 




Ay* (s) = -LI V ,\ xo^ 24s - 36m 
Bl+*~(s) = 0; = 2.2 x 10~ 3 . (2.21) 

The decay width from the amplitude ( |2.21| ) turns out to be very small (r(L 4 , tc + 7t~) = 



0.45 x 10~ 4 eV). In superposition to the amplitude corresponding to the decay width ( |2.20| ), 
its effect is totally negligible. 

The calculation of charged kaon loops from £ 4 is analogous to the pion loop calculation 
but invokes the terms in the Gasser-Leutwyler Lagrangian proportional to L4, L5, L§ and 
Liq. By interference of all terms, we gain a factor of 25 (T(L 4 , ti + tt~ + k + K~) = 1.9 x 10 -3 
eV), the complete £ 4 -contribution reduces the decay width ( |2.20|) to 



The values quoted for the £ -loops seem to match the Z A = 0.62 results of Ref. p2| . We 
come to the conclusion that the one-loop contributions from C 2 and C A mostly cancel each 
other. 

The next group of diagrams are the factorizable 2-loop diagrams derived from Fig. ^ (a) 
and (b). Again, we first consider the case where only pions are propagating in the loops. If a 
loop has no photons attached to it, there are also neutral pions allowed. Thus, the 6-meson 
vertex diagram of Fig. |3] (a) contains two charged pion loops plus the combination of one 
charged and one neutral loop. 

It is a straightforward procedure to determine the Feynman diagrams and derive the 
corresponding integral expressions. Solving these integrals is largely simplified by the fact 
that they effectively reduce to a product of two one-loop integrals. On the other hand, as the 
masses of external and loop particles are comparable, we need to solve the integrals exactly. 
This can be done by putting propagators with different loop momenta together by means 
of Feynman parameters and performing the four dimensional integration in the scheme of 



dimensional regularization O ne is left with an (at most) two-dimensional integration 



in the Feynman parameters space that is analytic due to the fact that the external photons 
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are on their mass shell. It is further simplified because the masses of the loop particles are 
all equal. The task is nevertheless a bit tedious because the one-loop integrals generate 
divergent terms proportional to lim„^ 4 if the number n of space-time dimensions is 
taken to 4 so that one is forced to expand all the expressions to linear order in e = 
Details of the calculation are discussed in the Appendix, subsection B. 

It is convenient to express everything as sums of products of the elementary integrals 
had, hub etc. (Eqs. (|4.2|), (|4.3|) , (|4.6|) ) and then do the multiplication of one-loop integrals 
and subtraction of singularities numerically. In order to be consistent with the scheme of 
chiral perturbation theory and the renormalization of the C 4 coupling constants, we have 
to subtract terms proportional to powers of R= — \ — 1— ln(47r) + 7, where 7 is the Euler 
constant. Also, the renormalization constant is to be chosen as \i — m v in order to be 
consistent with the numerical values quoted for the £ 4 low-energy constants. The finite 
end result has to be gauge invariant, i. e. proportional to g^ v \ — q^Qi- (As in the case of 
the one-loop amplitudes, there is no contribution B(s,t).) Gauge invariance can be shown 
analytically and, on the other hand, serves as a check to the calculation. 

It turns out that the various terms in the factorizable two pion-loop amplitude largely 
cancel each other. Our result for the decay width is T = 2.7 x 10~ 4 eV, about a factor of 3 
smaller than the one pion-loop amplitude from C 2 and well in agreement with the result of 



The effect is also very small in superposition 

r Le + o( P 4 )+L4+fact.^ ^ = _ 682e v. (2.23) 

Due to the crucial role of cancellations between various terms in the amplitude, it is 
problematic to predict contributions without actually performing the calculation. We may, 
however, observe the following: Replacement of a pion loop by a kaon loop lifts the G- 
parity suppression expressed by the factor (Am^) 2 /m 2 K but introduces a mass factor from 
the propagators which, in the low-energy limit, amounts to m 2 /m 2 K . A very rough estimate 
would therefore be that the amplitude from a kaon loop is larger than the corresponding pion 
loop amplitude by a factor of m 2 /(ArriK) 2 ~ 3, a factorizable two-loop diagram containing 
two kaon loops should be suppressed by m 4 r /(Am^) 2 m| : ~ 1/3. The values 3 and 1/3 
quoted are the relevant ratios because only the interference terms between loop and tree- 
level contributions affect the result for the decay width noticeably. Our estimate about 
correct for the C 2 one-loop contributions and seems to be confirmed by the two-loop results 
quoted in [12|. (Estimating the kaon loops from £ 4 is more problematic because more 



structures in the Lagrangian can generate kaon loops.) Applying it to the factorizable two- 
loop amplitude, we would conclude that the two-kaon loops should be totally negligible 
whereas the pion-kaon loop amplitude might change the result ( |2.23p by ±0.02 eV, a value 
that is markedly lower than the uncertainties introduced by the meson exchange model. 

The technically challenging part of the amplitude consists in evaluating the overlapping 
two-loop diagrams Fig. |7|. Again, we need the exact result and can not resort, e. g. to 
developments in the external momenta that have been successfully applied for QCD calcu- 
lations [SSI]. Feynman diagrams of the type required have been evaluated for the 77 — > 7r°7r° 
amplitude. As outlined in Ref. the expressions corresponding to Fig. |7| (a) and (b) can 



indeed be transformed into two-dimensional integrals by first integrating over one of the 
loop momenta and then transforming the result into a dispersion integral over a parameter- 
dependent box diagram that in turn can be represented as a dispersion integral. We note 
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that the graph [7| (a) leads to direct s — t channel box diagrams suggesting a fixed-i dispersion 
relation whereas the "master" -diagram [7| (b) yields crossed u — t channel box diagrams. In 
our case, the analyticity properties of the graphs are unfortunately more complex because 
the external 77 is unstable against decay and thus introduces anomalous thresholds |22| . Cor- 
respondingly, we will find that the overlapping amplitude has a complex part for arbitrary 
provided that > 3ra„-. 

Our solution avoids dispersion relations and is inspired by treatments such as where 



the graphs are transformed into expressions of the type I(m,n,Z) (see Eq. ( |4.12|) ). The 
hierarchy of I(m, n, Z) is interrelated by partial fractions, differentiations and partial inte- 



grations so that, in the case k = considered in Ref. ||23|| , all integrals can be derived from 
the basic formula for 1(2, 1,1) which itself can be written as an analytic expression. 

For the problem considered here, we have to extend the formalism by introducing a four 



momentum vector k in the expressions ( [4.12 ) which depends on the external momenta plus 
up to 2 additional Feynman parameters. This implies no major changes in the Feynman 
parametrization and evaluation of the four momentum integrals but markedly complicates 
the subsequent Feynman parameter integrations: in general, we obtain a four dimensional 
integral where, due to the 77 instability mentioned above, the integrand contains logarith- 
mic singularities within the parameter domain and is complex for arbitrary values of the 
kinematic variables. It is essential to work out the analytic structure of the integrand by per- 
forming at least two integrations analytically; the remaining two-dimensional integral can 
then be estimated numerically with reasonable accuracy. We will point out the nontrivial 
steps of the procedure in the appendix, sect. C. 

Because of the amount of algebra involved, it is essential to do the whole calculation as 
systematically as possible. As far as the derivation of integral expressions and the reduction 
and simplification of the numerators Z(li, l 2 , q±, q 2 , P,Xi,---, x m ) is concerned, we are helped 
by the fact that the 7737T- vertices always produce a constant that factors out. The basic 
integrals I(m,n, Z) can be checked independently and various relations between them help 
to exclude errors in the derivation. Finally, at least in the subthreshold region, all analytic 
one- or two-dimensional integral formulas can be easily checked on the computer. 

Unless the loop corrections discussed so far, the overlapping diagrams contribute to the 
amplitude part B(s,t) as well as to A(s). This implies that there is a contribution from 
these graphs even at s = where A(s) has no influence (see Eq. ( |2.7| ) and Fig. |§p. The 
analogy with the pion production process where nonanalytical contributions were found to 
be small in the energy range 2777^ < y/s < 400 MeV |U| can therefore not be used to argue 
that the effect of the overlapping diagrams on the //-decay rate is negligible fll2"] . Indeed, we 



find a decay width from the overlapping pion diagrams alone that is comparable to the one 
found for the factorizable diagrams r = 4.1 x 10~ 4 eV, this changes the decay width (|2.23 ) 
to 

r L»+0( P *)+I,*+fect.+ovBrL( 7? ^ = . 6 96eV. ( 2 .24) 

Again, we might estimate an uncertainty of ±0.04 eV due to the kaon loops. If - as in the 
case of the 0(p) 4 one-loops - the interference of kaon loops and pion loops is equal, the 
prediction of chiral perturbation theory to 0{p) & should be am 77-decay width of about 0.76 
eV. 
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As pointed out in Ref . j| , the anomalous part of the £ 4 Lagrangian generates one-loop 
terms that are 0(p) s in the momentum expansion but of a size comparable to the 0(p) 4 
loops because there is no G-parity or kaon mass suppression. Adding these contributions 
(Eq. (27)/(28) of Ref. @) to the amplitude of Eq. (ggg ), we get the final result 

rL 6 +0(p 4 )+L 4 +fact . +overL ^ ^ = q 765e y ( 2 _ 25 ) 

The coupling constants of the vector mesons used to determine the £ 6 -constants are 
subject to errors. Based on the uncertainties quoted in || for the u and p couplings, we 
calculate a relative error of the tree level amplitude ( [2.1 9|) of about 20 %. Together with 
the missing kaon loops, this implies an uncertainty to the 0(p) 6 amplitude of 

rL 6 +0(p4)+L 4 +fact . ^ ^ = 0.77 ± 16 eV. (2.26) 

This is in close agreement with the experimental value and the result of the quark box model 
calculation of [2^] but is markedly higher than the prediction of the meson exchange model 
used to determine the Lagrangian. 

Addition of the loop contributions to the "all order fit" amplitude of Table II does not 
alter the interference scheme described so far but would result in a total decay width of 

pVMD+loops^ _^ = 439 ± 9e y_ ( 2>27 ) 

The difference between Eqs. (|2.25|) and ( |2.27| ) gives an estimate of the uncertainty of the 
model induced by higher order terms in the momentum expansion. The various contributions 
to the decay width in both the 0(p) 6 and the all order fit scheme are collected in Table 3; 
the influence of the loops to the decay spectrum is shown in Fig. [8|. The plot shows again 
that contributions to B(s,t) affect the spectrum especially at low energy whereas the A(s)- 
contributions are cut off at s = 0. 



III. CONCLUSIONS 

We have considered the decay t] — > 7r°77 as an example for the application of the general 
chiral Lagrangian £ 6 . A total of 6 structures contribute to the process; the corresponding 
low-energy constants can be determined in a meson exchange model if the analogous process 
77 — > 7r°7r° is considered simultaneously. To be consistent at 0(p) 6 , the £ 6 tree-level ampli- 
tude must be complemented by one- and two-loop amplitudes. The result is in agreement 
with the experimental data for both the neutral pion production and the 77-decay process. 

We have given arguments for how to choose the relative signs of the coupling constants 
in the meson exchange amplitude. As this amplitude contains terms of 0(p) 8 and higher, 
there remains a sizeable difference between the predictions of the full meson exchange model 
and 0(p) 6 chiral perturbation theory. The total decay widths calculated in both models 
differ by a factor of 1.7 (F(VMD) = 0.45 eV, T(xPT) = 0.77 eV. Furthermore, unlike the 
meson exchange spectrum, the differential width in chiral perturbation theory is small at low 
invariant energy s and reaches a maximum at about 80 % of the maximum s kinematically 
allowed. An experiment measuring the decay spectrum would discriminate between both 
predictions and thus shed light on the details of the mechanism. 
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The loop contributions for the 77-decay are all individually small. However, as we know 
from the contribution of the (9(p) 4 -loops there can be sizeable interference effects with 
the leading tree-level amplitude. After a detailed calculation we could show that such 
interference effects occur (at the 5-10% level) in parts of the amplitude but that the complete 
two-loop result is suppressed by destructive interference, in particular between the £ 4 kaon- 
loops and the rest. The two kaon-loop amplitudes have not been calculated but rather 
estimated to be sufficiently small. Our results confirm the convergence of the chiral loop 
expansion. 

The relative size of the two-loop contributions, in particular the overlapping graphs 
displayed in Fig. [7] are not quite large enough to justify the amount of work required for 
their calculation. Nevertheless, we are confident that the strategy for solving the graphs 
outlined in the appendix may be useful for related calculations. In particular, it would 
be worthwhile to treat the production amplitude 77 — > 7r°7r° in SU(3)-chiral perturbation 
theory and compare our formalism with the alternative dispersion-theoretical approach. 

IV. APPENDIX: A RECIPE FOR SOLVING TWO LOOP INTEGRALS IN X PT 



In this section, we develop a method for solving the two-loop contributions for a two 
neutral meson/two photon process applicable for equal masses of the loop-mesons and real 
photons but arbitrary external mesons. (The general mass case would only introduce some 
modifications in the analytic integrations of the overlapping diagrams treated in subsection 
C.) The general diagrams contributing at 0(p 6 ) in the chiral expansion are shown in Fig. 2 
of Ref. ||13|| . For the //-decay, the diagrams involving £ 4 tree-level interactions are cancelled; 
the only diagrams involving £ 4 -terms are thus one-loop diagrams with a single 4-meson 
vertex. 

As C 2 can only generate vertices with an even number of mesons, the two-loop configura- 
tions we are left with are those displayed in Fig. [| where the squares denote £ 2 -interactions 
and photons are to be attached in all possible ways. If we restrict ourselves to internal 
pions, the diagrams contain each a vertex function Y vlcnv =const. (see Eq. ( |4.1| ) below). By 
a symmetric loop integration argument, one can show that this excludes the diagrams of 
the type (c) of Fig. f|. We are thus left with the factorizable diagrams (a) and (b) plus the 
overlapping diagrams of type (d). 

The vertex functions required can be calculated in a standard way by inserting the 
expansion U — 1 + — . . . in the Lagrangian C 2 and collecting all terms with a given 
number of meson and photon fields. The procedure implies the trace of a product of up to 
six 3 x 3-matrices and is therefore most conveniently done on the computer. Neglecting the 
kaon sector, one obtains the following vertex functions: 



A. Restrictions, Remaining Vertices and Diagrams 




(Am K ) 2 



14 



rw-^- = ^ (p\ + p 2 2 +p 2 + + p 2 ^- 3{p+ - P-f + 2m] 

2ie 8i 

— nW. V^ u — n tiv. 



^p2 iJ ' 7T+7T 7T + 7T~ 3^"^ 

z \/2 /a \ 2 ^ 2 .■> 



All other vertex functions vanish. The conventions used for these expressions, in particular 
the orientation of the four momenta, are displayed in Fig. |5|. Statistical factors for the case 
that there are two or more identical particles have been included in Eqs. (}4.1|) . For vertices 
involving two 7r + 7r~ pairs, the momenta p + , p_ have to be assigned symmetrically. Note that 
there are no photon vertices with neutral mesons only and that all the vertices involving 
a single ^-meson are generated by the mass term of the Lagrangian C 2 and are constant. 
In analogy to the one-loop amplitude |3J, the total internal pion two-loop amplitude is 
G-parity suppressed. 



B. Factorizable Diagrams 



The Diagrams (a) and (b) of Fig. f| contain no propagators containing both loop momenta 
and therefore effectively factorize into a product of two one-loop expressions. According 
to the scheme of dimensional regularization [^T|, one-loop integrals contain singularities 

linear in - = — where n is the continuous dimension parameter. In order to retrieve all 

6 1 2 

nonvanishing terms as n — > 4, we are forced to develop the one-loop expressions to O(e). 
The divergence of the resulting two-loop amplitude is at most quadratic in -. 

For clarity, we display the explicit expansions of the simplest elementary one-loop dia- 
grams. A tadpole loop (see Fig. || (a)) generates p3| : 



Had 



dH 1 



im 2 T(2 



z 2 



(2tt) 4 [/] (4tt) 2 1 



H 

2 



77T 



im 



(47T) 



+ 



m 



1-7-W -) 



m 



2 1 2 
m 1 2/ m - 

+ 2 n ^4^2' 



+ 0( 




(4.2) 



where we have denoted the propagator as [/] = I 2 — m 2 + ie, expanded the T-function 



r(e) = - - 7 + eS + 0{e 2 
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with the Euler constant 7 = 0.577215665 and introduced the renormalization constant fi. 
The small imaginary factor in the propagator will not be written explicitly in the following 
but should be kept in mind in order to retrieve the imaginary parts of the loop integrals. 
Also, in Eq. (|4.2j) as well as in the following, the limits n — > 4 and equivalently e — > are 
always understood. 

A symmetric bubble-diagram (Figure ^ (b)) with photon momenta q\ + q2 where s = 
{Qi + Q2) 2 yields after Feynman parametrization 
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1 bub 



(2 7r )4[/ + 2i±22][/_2i±2a] (4tt) : 



9 n 

Z 2 



(1 — x) — m 2 



(4vr 
+e 



+ 



7 -ln(— -G(s)) 



"47T/i 



5 + ^^4V ) + 2 ln( 



+ 7 G(s) + In 2 



m 



rG(s) 



where the loop integrals 

G(s) = [ dxln(l - -4?ar(l -x))) G(0) = 0) 
Jo m z 

G 2 (s) = [ dxln 2 (l- -i-x(l-x)); G 2 (0) = 
Jo m 2 

are analytic functions cut along the real axis as s > Am 2 . The integrals G(s) and 

1,1 dx 



+ 0{e 2 ) . (4.3) 



F s) 



X 



ln(l 



r x(l-x)); F(0) = 



(4.4) 



(4.5) 



can be easily reduced to analytic expressions |I|,|13|; moreover, loop integrals with a single 
logarithm but additional factors x, x 2 etc. can be expressed in terms of G(s) and F(s). The 
loop functions of the form G 2 (s), F 2 (s) etc. can be expressed in terms of polylogarithms 
f2~4fl or be solved numerically. 

The scalar vertex diagram (scalar meaning numerator equal to 1, see Fig. ^| (c)) is finite 
and written as 



dH 



(27r) 4 [/][/- gi ][/ + g 2 ] (4vr) 2 Jo Jo 



dx\ 
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dx 2 T(3 ) 



(AtT/J, 



2\2- 



(4tt) 2 Jo 
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CLXo ~ — 



2 (sx\X2 — m 2 ) 3 2 



3X1X2 



sxix 2 — wi z 



(4.6) 



The scalar integrals considered so far can be easily generalized to the case where loop 
momenta appear in the numerator of the integrand. For the case of a linear momentum 
in the numerator, the formula for dimensional regularization are applicable. Integrals with 
quadratic and higher order loop momenta in the numerator are reduced to simpler integrals 
with formula such as 



dH A\n v - V^ 2 



AWk" - ^g^k 2 



(2tt) 4 (Z 2 + 2kl-m 2 ) a 16ttS (-k 2 - m 2 ) 



(4.7) 



The vertex and box diagrams appearing in the amplitude generate further loop functions 
of the types 



dx\ 



1— XI 



dxo 



x™xV^ 



rl rl—Xi 

/ dx\ 
Jo Jo 



dx? 



SX1X2 — m 



1 Jj 2 
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dx\ 



l-x 



1 x™X2 ln(l — ^2-XiX 2 ) 



dx 



(sxix 2 — m 2 ) 2 
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with powers < m,n < 3. The loop functions are symmetric in m and n. In all cases, 
the ^-integration is easy to perform analytically; partial integration and simple algebraic 
manipulations allow to express everything in terms of F(s), G(s), F 2 (s) (defined as F(s) 
but with a quadratic logarithm in the integrand) and G 2 (s). Due to the similarity of the 
loop functions, we find that the factorizable amplitude has a unique s-channel threshold 
s = 4m 2 . 

Let us illustrate the method in a simple example, the 6-meson vertex diagram of Fig. |4] 
(a). In order to get a nonvanishing amplitude, both photon lines have to be attached to one 
charged meson loop; the second loop can be charged or neutral. From the vertex expressions 
in Eq. QL1|), one finds that the amplitude can be written as 



3tV2e 2 {Am K ) 2 r d% d% {2k - q x Y{2h + q 2 ) u ~ 9^(11 



2\ 



M " "2 9^ 7 (2^(2^ [/i] [/! + *] + 

(4.8) 

This is gauge invariant because of 

rf 4 /i qi, [(2h ~ qi)»{2h + <li) v ~ 9^{l\ ~ ™ 2 )} 



qi,M^ 



(2tt) 4 [h}[h-qi}[h + q2\ 

d% 2l\ r d% 2l\ 



( 2vr )4 [ h - ffl±E] [ Zl + 2i±£a] J ( 27r )4 [ Zl - f ] [I, + f 



0. (4.9) 



With the expressions ( |4.2|) , (|4.3| ) and (|4.6|) for the elementary diagrams, the amplitude ( [4.8|) 

reads 

— -^f^( 9 i-«)^x 

i f . ^ ^ [i + ; (-7 - Mg) - Hi - **,*,))] _ 



(4ir) 2 Jo Jo sxix 2 — m? 

The loop integrals are readily solved and yield 



1 f 1 -^ , XiX 2 l/l m 2 
dxi / dx 2 5- = - - H 

JO SX1X2 — 777 s \2 s 



X j f l ~ xi j ^1^2 ln(l — ^a;ix 2 ) 1 / 1 m 2 J71/ N i 1 nll s i m 2 ^ 



^1/ da*-^^ = - + -G(a) + — F 2 (s) 

Jo sxix 2 — m 2 s \ 2 s 2 2s J 

(4.11) 

The 6-meson diagrams contribute to the form factor A(s,t) of Eq. ( |2.6j ) only. The same is 
true for all factorizable diagrams. 



C. Overlapping Diagrams 

The overlapping diagrams derived from Fig. ^ (d) translate into sums of integrals of the 
type 
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rl rl—xl— ...—Xn-l rl—Xn+l — — —X m -1 

I(m, n, Z) = / dxi ■ ■ ■ dx n / dx n+ \ ■ ■ ■ dx m 

Jo Jo Jo Jo 

[ d h d l 2 Z(li,l 2 ,q±,q2, P,Xi, ■ ■ ■ ,x m ) (a-\o\ 
X J (2tt)* (2vr)4 W\QF\h ~k- Hx lt ~ x m )] ' [ ' } 

where 1 < n, m — n < 3; m < 4. Here, we have used Feynman parameters in order 
to put together propagators with equal loop momenta and obtain more symmetric expres- 
sions. The numerators Z can be largely simplified by factoring out /-dependent terms; the 
remainder is at most an expression cubic in the loop momenta \H v l p . With further Feynman 
parametrization, all the integrals I(m,n, Z) can be cast in a form that makes the analytic 
structure of the loop integrals explicit and leaves a two-dimensional numerical integration 
over a well-behaved, complex function. We demonstrate the method for the scalar integrals 
7(2,1,1), 1(2,2,1) and 7(3,1,1). 

The integral 1(2,1,1) corresponds to the scalar part of the diagram Fig. |^ (c). The 
analytic case where k = has been considered in ||23|| . In extension to the method outlined 
there, we apply subsequent feynman parametrization and integration over the loop momenta 
and write 

- 1 , r d% d% 1 



1(2, 1, 1) = J Q dxJ 



(2tt) 4 (2^ [h] 2 [l 2 ][l 2 - h - k(x)] 
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{ k2 2fc2 (! -y) + k\i- yf) r(5 - i 



(k 2 y(l -y)~ m 2 y - m ^_ z y ) ) ) 
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+2 Eii^L____^ (l .i 3) 

2 (k 2 y{\ -y)- m 2 y - "ffj^ ' 

After expanding (1 — y) 2_ ^| in powers of e = 2 — |, the second term in this expression 
can be partially integrated over y. One obtains a sum of terms where the z-integration plus 
some trivial x- and y-integrations can be done analytically. The final result can be cast in 
the form 

-1 f 1 if ( m 2 \\ 

1(2, 1, 1) = — — — + - 1 - 2 T - 2 In 



(4tt) 4 [2e 2 2e \ 1 \4ir{i 2 J J 

+ - + --7 + 7 2 -(l-2 7 )lnf^)+21n 2 

nl dyG(a)- [ 1 ^(G(a) + Hy) + 2)\; 
) Jo y J 

k 2 (x)y(l -y) -m 2 (l-y) 
m 2 y 



+ 



(4.14) 



with the loop function G from Eq. ( |4.4j) . We note that for the 77-decay diagram Fig. |7| (c), 
the external momentum variable can be written as k = P — q±x . Owing to the properties 
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of G and the ratio of the 77- and 7r-masses, the integrand in Eq. |4.14] is a complex function 
for arbitrary s > 0. The integral is finite with a logarithmic end point singularity in the real 
part of the integrand as y — > 0. 

Feynman integrals with a numerator depending on the loop momenta can be solved 
analogously. One obtains 



1(2,1,^) EE J^dxJ 



d% d% It 



(4.y r { 47 - m - ^ (4$ ) + r ^ - ^ 

J(2,l,A/n = / o ^7^)4(2^)4 [Zl ]2 [Z2 ] [Z2 _ Zl _ A . (a . ) ] 
= — L dx(4fc(x)"fc(x)" - -*(*) V") 

(47TJ 4 JO n 

(4.15) 

The last equation shows how to reduce integrals with higher powers of loop momenta in 
the numerator to the integral class 1(1, 1, Z). The latter integral can in turn be related to 
1(2, 1, Z) by means of the 'partial' operation JO],^. 

The 'master' diagram Fig. [F| (b) leads to a scalar integral 7(2,2,1) with k(xi,x 2 ) = 
P — qiXi — q 2 x 2 - The integral is finite so that we don't have to worry about the regularization 
procedure: 
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= (ST X #1 k dm k dZl k dZ2 Pz^ + Pz^-zJ-m* 

=w?^r w ( f i^f) - f )) • (4 - ie) 

where r 2 , s 2 and fh 2 depend on the Feynman parameters y±, y 2 and the invariant kinematic 
variables 



r =sy 1 y 2 -{m -u)y 1 (l-y 1 ), 
s 2 ee -(m 2 - t)yty 2 - m 2 (l - yi) 2 , 
m 2 ee m 2 yi. 



The analytic integration is largely simplified by the fact that both photons are real 
(Qi = Q2 = 0) an( i the masses in the propagators are identical. The general mass case 
can still be treated in the same way, but with polylogarithms occurring from the analytical 
^-integration. 

The integral 1(2, 2, 1) is linked to 1(2, 1, 1) by 
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7(2,2,1)= lim 
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/ dx\ I dx 2 
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(4.17) 



Relations such as Eq. ( [4.171 ) can De used to check the calculation. 

For an accurate numerical estimate of two-loop Feynman integrals, it is crucial to evalu- 
ate the first complex integration analytically. Due to the x-parametrization of loop momenta 
and the presence of /-dependent terms in the numerator Z, the master diagram contains 
integrals of the form 



I M (l,m,n) 



dz\ 



dz-2 



z\ z. 
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1 - z 1 ) n aziz 2 + 62:1(1 - z\. 



0</,m,n<3. (4.18) 



The case m > leads to cancelling singular terms after the first integration and requires 
formula such as 

H f {l + f z (in (l " " z)) - In (l - \z{l - z)))} 



that are conveniently derived by formal partial integration of the logarithmic terms and 
subsequent rearrangement of the integrands keeping in mind the infinitesimal complex in- 
crement so that, e. g. 

2 (G(a) + 2). 



lo az 2 — az + 1 4 — a 

It turns out that all integrals I M (l,m,n) can be expressed in terms of the loop functions 
F and G. As in Eq. ( f4.14|) , we are left with a two-dimensional integral over a complex 
function that can be easily done numerically. 

The diagram Fig. [7] (a) yields k = P — q\(l — x-y) — q 2 x 2 and the modified integral 



7(3,1,1) 



dx\ 
-1 rl 



1— x\ 



dx 



d% d% 



1 



(47T) 



dx\ 



2 

1— Xl 



(2vr) 4 (2vr) 4 {I 2 + sx lX2 - m 2 f[l 2 ] [l 2 - h - k{x u x 2 )} 

,2\2 \ 2 ~~ 







dx 2 
k 2 y - 







z(l-z) 



e \z(l-z) 
r(6-n) 



dyy 2 (l-yy * 



(k 2 y(l -y) + (sxix 2 - m 2 )y - m ^_^ 



n 
+ 2 



r(5-n) 



(k 2 y{l - y) 



[sx\x 2 — m A )y — 



m 2 (l-y)\ 5 - n 



z(l-z) 



(4.19) 



At this point, it is convenient to write the denominators as 



k 2 y(l — y) + {sxix 2 — m 2 )y — m ^ ^ 



z(l-z) 



1 i- d 

y(t-l) llmM2 -™ 2 9M 2 



k 2 y(l -y) + (sx lX2 - M 2 )y - 



-t+i 
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partially integrate the second term as in Eq. ( |4.14| ), expand and rearrange all terms and 
finally reperform the partial derivation. The result can be represented as 

— 1 f 1 rl-Xt ( rl rl dy 

1(3, 1, 1) = / dx\ / dx 2 I dz 



(47r) 4 jo Jo [Jo Jo 1 — y 

£ _ 2 i__ 

k 2 y(l -y) + {sx x x 2 - m 2 )y - s Xl x 2 - m? 

+ f d M^m>- \ r(5-n) ) 
Jo e \z[l — z) ) \sx\X2 — m*) n \ 

Note that 1(3, 1, 1) contains simple s-channel threshold terms in sx\X2 — m 2 as well as 
an anomalous threshold term. The cancellation of the parameter singularities in the first 
two contributions to Eq. ( |4.20| ) can be used as a test for the analytic integration. The 
x-integrations in the s-threshold terms can be performed according to subsection B. In the 
first term of Eq. (|4.20| ), the x%- and ^-integrations can be done analytically yielding a form 
analogous to Eq. ( |4.16| ). As in the case of the master diagram (see Eq. ( J4.18[ )), the diagram 
|7| (a) involves a whole set of analytic parameter integrals. 

In analogy to the integrals of the class 1(2, 1, Z), singularities may cancel in the presence 
of /-dependent numerators. For instance, we obtain 

1 , f 1 -^ j f d% d% 11 



7(3, 1, Zi) = / dx\ \ dX2 

Jo Jo J 



(21)" (2tt) 4 ((? + sift - m 2 ) 3 [i 2 ][i 2 -h- k(x)] 
Jo Jo Jo Jo 



~ (4vr) 4 Jo Jo Jo Jo y k 2 y(l - y) + [sx 1 x 2 - m 2 )y - 

Jo Jo J 



(2tt) 4 (2tt) 4 (I 2 + SXl X2 - m 2 )3[/ 2 ][/ 2 - I, - k(x)} 

— Cdx C' Xl dx CdzCd y 2 (l-y)(k^- g ^k 2 ) 

(47r) 4 io Xl Jo X2 Jo Z Jo V k 2 y(l - y) + (SX.X2 - m 2 )y - 
This closes the list of basic integrals required for the type of processes we are considering. 
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TABLES 



Meson m [MeV] G w [GeV -2 ] G v [GeV~ 2 ] L^eV] 



p 


768 


0.084 


0.17 


0.11 




782 


0.49 


0.098 


0.28 




1019 


0.018 


0.0089 


0.29 


bi 


1232 


0.39 


0.18 


0.32 


a 


983 





-0.018/-0.022 


0.35 




1318 





0.216 


0.37 


h 


974 


-0.007/-0.009 







h 


1275 


1.28 








TABLE I. Meson parameters used for the meson exchange amplitudes of Eqs. ( |2.15| ) - ( [2, 18 ). 
For the scalar mesons, both coupling constants G Sd and G Sm are given. The last column contains 
the (cumulative) total ??-decay width. 



Model 


di 




d 2 


d 5 


d 3 


d.6 


0(p 6 ) 
All Ordei 


-1.60 
0.108 


0.654 
1.64 


3.80 
0.0486 


-1.93 
-4.13 


38.9 
40.3 


-16.1 
-16.7 


TABLE II. £ 6 low-energy constants determined in two different fit schemes (see text). All 
constants are in units of 10 -3 GeV~ 2 . The relation between d\,. . . ,d§ and the low energy constants 
of Ref. |7| are given in Eq. (2.12]) 


Model 


£ 6 0(p) 4 


£ 4 


fact. 


overl. 


anom. 


(K-loops) 


0(p 6 ) 
All Order 


0.652 0.733 
0.371 0.411 


0.673 
0.381 


0.682 
0.383 


0.696 
0.386 


0.765 
0.439 


(0.83) 
(0.46) 



TABLE III. Cumulative contributions of the loop contributions to the total r^-decay width. 
The £ 6 -fit parameters are those given in Table II, the loop corrections are denoted as in the text. 
The last column contains the estimate of the missing 0(p) 6 kaon loops. 
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FIGURES 



FIG. 1. The 77 — > 7r°7r° spectrum as measured by the Crystal Ball Collaboration fl20f . The 
curves are results of chiral perturbation theory complemented by a 0(p G ) meson exchange ampli- 
tude: Pure chiral one- and two-loop contributions in the SU(2)-model (Eqs. (7.2) and (7.7) of Ref. 
p~3[| ) (dashed), loop amplitude plus full meson exchange model Eqs. ( |2.15| ) - ( 2.18 ) (dotted) and 



loop amplitude with meson exchange model but without axial mesons Eq. ( [2.16 ). 



FIG. 2. Predictions for the r/ — ► 7r°77 spectrum. The predicion of the full meson exchange 
model (solid) and its best fit according to the amplitude Eq. ( 2.13|) (dotted) yields a total decay 



width of T to t = 0.37 eV. For comparison, the predictions of the 0{p ) fit (dashed, T to t = 0.65 eV) 



and the 2-parameter fit according to ||], Eq. (LI) (short dashed, T to t = 0.18 eV) are also shown. 
Loop corrections are not taken into account. 



FIG. 3. 77 — ► 7r°7r° cross section from [2C] and predictions from chiral perturbation theory. ^. 
The results have been obtained with the chiral one- and two-loop contributions (Eqs. (7.2) and 
(7.7) of Ref. |L|) complemented with the 0(p e ) -contributions as in Fig. |. Axial vector mesons 
are not taken into account. 



FIG. 4. The two types of factorizable two-loop diagrams contributing to the r\ — > 7r°77 am- 
plitude and considered in this paper. 



FIG. 5. Vertices and vertex functions r^ +?r _, r^ 07r07r+7r _ and T 7T + n - 7T o n o w o ri , as used in Eq. (4.1). 



FIG. 6. Basic one-loop diagrams calculated in Appendix, subsection B: Tadpole (a), bubble 
(b) and vertex (c) diagram. 



FIG. 7. The four types of overlapping diagrams occuring in the process rj — > 7r°77. 



FIG. 8. Effect of the loop contributions on the 77-decay spectrum. The £ 6 tree-level results 
according to Table III, first column are solid (compare Fig. Q), the dashed lines include the 0(p) 4 
one-loop result (Table III, second column), the dotted lines correspond to the full result (Table 
III, second last column). 
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